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Exotic stable massive particles (SMP) are proposed in a number of scenarios of physics 
beyond the Standard Model. It is important that LHC experiments are able to detect 
hadronic SMP with masses around the TeV scale. To do this, an understanding of 
the interactions of SMPs in matter is required. In this paper a Regge-based model of 
_R-hadron scattering is extended and implemented in Geant-4. 
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1 Introduction 

The observation of exotic stabl^l massive particles (SMPs) with electric, magnetic or colour 
charge (or combinations thereof) would be of fundamental significance. Searches are there- 
fore routinely made for SMPs bound in matter, and for SMPs at cosmic ray and collider 
experiments [U [21 [3l S] . Furthermore, SMPs with masses at the TeV scale are anticipated 
in a number of beyond-the-Standard-Model scenarios, such as supersymmetry and universal 
extra dimensions . SMP searches are thus a key component of the early data exploitation 
program of the LHC experiments [S]. To facilitate such searches, models are needed of the 
interactions of SMP as they pass through matter, since such interactions can give rise to a 
range of different event topologies j4|. Although robust phenomcnological approaches are 
available to predict the scattering of electromagnetically charged SMPs, it is not clear how 
the interactions of hadronic SMPs (hereafter termed i?-hadront0) should be treated owing 
to uncertainties associated with the strong scattering processes and the hadronic mass hi- 
erarchy. This paper describes an implementation within Geant-4 [5] and extension of a 
Regge-based model [7] for hadronic scattering. The paper is a summary of work contained 
in Ref. [8]. 

To aid the development of search strategies, it is important that experiments have access 
to Monte Carlo models which span as fully as possible the range of conceivable signatures 
which could be associated with i?-hadron production. At present, one model is available 
within Geant-4 [3]. This model, hereafter termed the generic model, is based on a black 
disk approximation [1^ and employs pragmatic assumptions regarding possible scattering 
processes and stable i?-hadron species. This work concerns an approach, hereafter termed 
the Regge model, which uses triple regge formalism to predict the scattering of i?-hadrons 
formed from exotic colour triplet squarks (stop- like (t) and sbottom-like (6)) and electrically 
neutral gluino-like {g) colour octet particles. Furthermore, the model uses a different set 
of well motivated assumptions regarding i?-hadron mass hierarchies than those employed in 
the generic model. 



''The term stable is taken to mean that the particle will not decay during its traversal of a detector. 

''The term ij-hadron has its origin in the i?-parity quantum number in supersymmetry theories and the 
work in this paper will be presented in the context of supersymmetric particles. However, the results are 
generally applicable to stable heavy exotic hadrons. 
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2 Heavy hadron mass hierarchy 



Although there exist no complete calculations of expected i?-hadron mass spectra, estimates 
have been made using a variety of approaches such as the bag model [TTl [121 US] and lattice 
QCD |14| . These predictions, together with with measurements of heavy hadron masses, 
make it possible to approximately determine those features of the mass hierarchies which are 
most relevant for the modelling of i?-hadron scattering in material. Of particular interest 
are the masses of the lowest lying states, to which higher mass particles would be expected 
to decay before being able to interact hadronically. 

Given the observed spin and flavour independence of interactions of heavy quarks with 
light quarks |15] the lightest meson and baryon states can be inferred [16| from the mea- 
sured mass spectra of charm and bottom hadrons. Thus, charged and neutral squark-bascd 
mesons would be approximately mass degenerate and the lightest baryon state to which the 
other baryons would decay would be the scalar singlet qud i.e. the exotic equivalent of the 
Ac, A;, baryons. This picture is supported by calculations of hypcrfinc splittings of i?-hadron 
masses (TBI [TU] . 

The basic properties of giuino i?-hadrons may not be similarly inferred from Standard 
Model hadronic mass spectra and a greater reliance on phenomenological approaches is re- 
quired to estimate the stable states. Possible hadrons include giuino i?-mesons 
{guiP,gudP,gd(f) and gluino-gluon states (gg). Bag model [11] and lattice QCD [M] cal- 
culations estimate that the masses of lowest lying states of each of these hadrons differ by 
less than a pion mass. Detailed calculations of the lowest lying baryon states [121 US] using 
the bag model predict that the lightest state is the singlet guds. Consequently, higher mass 
states would decay weakly into the singlet state over a time scale of ~ 10~^°s [T5] . 

For the work presented here, we use the squark and giuino i?-hadron mass hierarchy 
assumptions given above and use a simplified mass spectra which comprises only those R- 
hadrons deemed to be stable. This leads to i?-hadron event topologies which would be 
challenging for the one of the most common collider search methods, which looks for a slow 
muon-likc object [3 [T71 [181 [191 [20] . This is due to the i?-baryon production processes which 
take place as i?-hadron pass through matter. Such processes, which, as discussed in the next 
section, imply that a sbottom or giuino i?-hadron, irrespective of its state as it entered a 
calorimeter, is likely to leave as an uncharged object and therefore escape detection in an 
outer muon chamber. 

2.1 Regge Model 

Since the central picture of a low energy light quark system interacting with a stationary 
nucleon [lOj . i?-hadron scattering can be treated with the phenomenology used to describe 
low energy hadron-hadron scattering data [71 [HI [^ , as is done in the Regge model of which 
more details can be found in Ref. [7] . The Regge model was originally developed to describe 
squark i?-hadron scattering though has been extended here to also treat giuino i?-hadrons. 
This model assumes the stable states described in the previous section. 

Using parameters fitted to low energy hadron-hadron data, the Regge model makes 
predictions for i?-hadron scattering cross sections, together with energy loss calculations 
based on the triple regge formalism. Fig. [T] (left) shows the the model predictions of the 
scattering cross sections of a squark-based i?-hadron off a stationary nucleon within a nucleus 
comprising equal numbers of neutrons and protons. The cross section is shown for different 
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Figure 1: Cross sections for stop-based and gluino-based i?-hadrons. The predictions are 
shown for the Regge model (sohd hnes) and the generic model (dashed lines). 



types of squark-based i?-hadrons as a function of the Lorentz factor 7. As can be seen, 
there is a large cross section for antibaryon (qud) interactions which is due to a dominant 
annihilation process with a nucleon in the target. 

It is also seen that, at lower values of 7 (7^10), the scattering cross section of squark- 
based i?-hadrons containing a light valence antiquark {qvP, qd) is larger than for antisquark- 
based i?-mesons, which arises from the presence of reggeon exchange processes which are 
only permitted for i?-hadrons containing a light anti-quark. Owing to the presence of an 
additional light quark the scattering cross section for i?-baryons (qud) is twice as large as 
for the mesons with light quarks. 

Upon an interaction, the probability that a i?-mcson to i?-meson process gives rise to 
charge exchange is 50%. Similarly, the probability of an interaction in which a i?-meson 
becomes a baryon is 10%. Once a i?-hadron becomes a baryon it stays in this state. Another 
process which must be taken into account is the oscillation of a neutral mesonic squark-based 
i?-hadron into its antiparticle. Since the conversion rate would be model dependent we allow 
two possibilities here: zero mixing, in which no oscillations take place, and a maximal mixing 
scenario in which there is a 50% probability that any neutral mesonic squark-based i?-hadron 
which was produced would automatically be converted to its anti-particle. 

In Fig. [1] (right) the predicted cross sections for gluino i?-liadron species are shown. The 
baryon cross section is 50% greater than that of stop i?-hadrons owing to the extra light 
quark in the gluino i?-hadron. Also, since the gluino meson contains a light quark and 
anti-quark, then both regge and pomeron exchange processes are available. Again, a baryon 
number transfer probability of 10% is used. A probability of | is assigned to i?-meson to 
i?-meson processes involving charge exchanges of 0, ±e or ±2e. 
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Figure 2: Number of hadronic interactions (left) and total energy loss (right) for i?-hadrons 
of 300 GeV mass sent through Im of iron. The top (bottom) plots show predictions for stop 
(gluino) _R-hadrons which entered the iron in the state td (gud). 



2.2 Implementation of Regge Model in Geant-4 

In the dynamical picture of i?-hadron scattering used by Geant, the light quark system 
is decoupled from the heavy spectator parton before interacting with a nucleon according 
to Geant's parameterised model for light hadrons. In this way secondary particles and 
so-called black tracks were generated. Following the interaction, the light quark system is 
recombined with the heavy parton. 

To implement the Regge model in Geant-4, the software architecture used for the generic 
model was adapted. The scattering cross sections and available baryon states were altered 
and the relative rates of the available processes were adapted to match the prescription in 
the above section. Using single particle simulations, it was verified that the main features 
(energy loss and charge exchange) of the Regge model were well reproduced in the Geant-4 
implementation. 



3 Interactions and Energy Loss 

To estimate the effects of scattering on i?-hadrons at the LHC simulations were made of the 
passage through iron of i?-hadrons of 300 GeV mass. The kinematic distributions of the 
i?-hadrons were given by the Pythia generator which simulated the direct pair production 
of tt, bb and gg in proton-proton collisions at 14 TeV centre-of-mass energy. 

The total energy loss for i?-hadrons passing through Im of iron is shown in Fig. [21 Stop 
i?-hadrons which start in the state td typically lose more energy than gluino i?-liadrons (gud) 
which start with the same charge. This occurs due to the rapid conversion of mesons to 
baryons. Since the gluino (stop) baryon is predicted to be of zero (positive) charge, it loses 
far less energy. Also shown is the number of hadronic interactions. Owing to the additional 
light quarks in the gluino i?-hadrons more scattering occurs in this case than for the stop 
i?-hadrons. A no-mixing scenario is assumed here. 
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4 Summary 



An implementation and extension of a Regge-based approach to squark and gluino-based 
i?-hadron scattering has been implemented in Geant-4. The work allows the modelling of 
the energy loss and the transformation of i?-hadrons as they propagate through material. 
Slides containing the presented work can be found in |23| 
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